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ABSTRACT 


Cross-shore wave transformation, nearshore currents, and morphology model 
predictions using Delft3D are compared with comprehensive observations acquired on a 
barred beach at Duck, North Carolina over a wide range of conditions. The Delft3D 2- 
DH model utilizes shallow water equations to phase resolve the mean and infragravity 
motions in combination with an advection diffusion equation for the sediment transport. 
Model coefficients and the effect of small changes in the wave incidence angle were 
examined for model sensitivity. The wave transformation model is dependent on the 
breaking parameter y, which determines organized wave energy dissipation. Y was 
found to increase as a function of offshore Hms. However, this is robust and a model 
skill of .89 was obtained using a constant y = .425. The manning number n affects the 
current bed shear stress and determines the model current magnitude having an optimal 
value of n = 0.02. The model is not overly sensitive to the value of n. The asymmetry 


coefficient a, determines the amount of onshore sediment movement. The rip channel 
created by mean currents on a short time scale is not affected by O,, whereas bar 
evolution requires more time to develop allowing o,, to affect morphology. A values of 


Q,, = .25 gave the best results. Further research is needed to calibrate this parameter. 


Small changes in wave angle can cause significant errors for currents when complex 
bathymetry is present and the waves are near shore normal. Overall the model is robust 


with sensitivity to small changes in near normal wave angles over complex bathymetry. 
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1. INTRODUCTION 


A. MOTIVATION 

Delft3D is being considered for transition as the new Navy surf model. Presently, 
the United States Navy uses the Navy Standard Surf Model, which is limited by 
simplified hydrodynamic equations (1-D) with the assumption of alongshore, uniform 
bathymetry. The 1D assumption simplifies the equations, but neglects the effects of 
alongshore variations on near shore dynamics. 2-D modeling is required for increased 
accuracy of morphology predictions in the nearshore. For example, 2-D model 
computations can reproduce the development of alongshore bathymetry of shoals cut by 
rip-channels (Reniers et al., 2002). Many areas such as coastal management and 
engineering rely on predictions made by 2-D models and would benefit from their 
improvement. Naval applications include amphibious, special operations, mine, and 
mine counter measures and would benefit from the increased accuracy of the 2-D model. 

The morphologic model that is examined is a research version of Delft3D. This 
2-D model utilizes the nonlinear shallow water equations to phase resolve the mean and 
infragravity motions in combination with an advection diffusion equation for the 
sediment transport (Reniers et al., 2002). Once Deflt3D free parameters are calibrated 
and the model skill is established, the model can help improve the knowledge of near 
shore dynamics. Beach morphology is intricate and must be studied with realistic data to 
improve the understanding of the dynamics. Roelvink and Broker (1993) state the need 
for detailed comparisons and verifications of sediment transport models. Initially the 
standard Delft3D model results were compared with the L1P11D experiment conducted 
in the Delta flume in 1993. Both erosive and accretive 1-D experiments were 
accomplished (Arcilla et al., 1994). The model was calibrated with the flume 
measurements and the results showed that hydrodynamic predictions performed better 
than morphodynamic predictions. Arcillia et al. (1994) state the need for additional 
modeling efforts, particularly for sediment transport modeling. There has been a lack of 
model verification using field data. The research version of the hydrodynamic model has 
been qualitatively compared with more complex 2-D situations in the RDEX field 


experiment in Palm Beach Australia (Reniers et al., 2001) and the RIPEX field 
1 


experiment in Monterey Bay, USA (Reniers et al., 2002). The morphodynamic model is 
assessed here using field data acquired during the Duck94 experiment at Duck, North 
Carolina, where an amphibious vehicle recorded almost daily variability of near shore 


morphology. 


B. OBJECTIVES 

Beach morphology is important in determining near shore hydrodynamics. 
Waves and currents vary as the near shore topography varies. Bars affect cross-shore 
wave transformation and alongshore bathymetry variations produce alongshore pressure 
gradients forcing longshore and cross-shore currents. The objective of this thesis is to 
assess the Delft3D model morphodynamics by comparing model results with Duck’94 
field data. 


Il. MODELING BAR MOTION 


A. MECHANISMS OF BAR FORMATION 

At least three processes have been identified that affect sandbar dynamics 
including wave asymmetry, undertow, and bound long waves. The break point 
hypothesis describes a mechanism of how sandbars are built. Outside the surf zone, there 
is an onshore movement of sediment by asymmetric wave motions. While inside the surf 
zone where waves break, the asymmetry decreases, and the set-up gradient drives an 
undertow, which transports the sand, stirred up by the wave motion, offshore. The wave 
asymmetry and undertow result in a convergence of sediment where the bar is located. In 
addition, wave group forced long waves outside the surf zone result in a seaward directed 
transport because of the concurrence of the bound-wave trough and the enhanced stirring 
during higher waves in a group (Roelvink and Stive, 1989). 

Wave asymmetry, a mechanism associated with bar generation, moves sand 
onshore. Owing to the inherent nonlinearities of shallow water waves, the crests of 
harmonic components become phase-locked resulting in an asymmetry of wave 
velocities. This asymmetry of the waves results in stronger onshore flow under the crest 


than offshore flow under the trough. Sediment transport, s, is a function of the velocity 
associated with short waves s ~ (u,?) where wu, is the intrawave short wave velocity. 


The result is sediment is moved onshore owing to the asymmetry of the wave. 

This onshore transport of sediments is countered by an offshore flux of sediment 
carried by the undertow. Waves produce an onshore mass flux due to wave drift and surf 
rollers confined primarily to the crest/trough region of the waves. This shore directed 
mass flux is compensated for by an offshore mean current below, the undertow. Wave 
breaking induced shear stresses are important in determining the vertical distribution of 
the mean velocity for undertow in the surf zone. The sediment stirred up by the wave 


motion is then carried off shore by this mean current. 


In addition to wave asymmetry and undertow Roelvink and Stive (1989) include 
wave group-induced long wave flow in the bar generation. Sediment is stirred by short 
waves on a wave group time-scale and is then subsequently transported by the 


infragravity orbital motion s ~ (usu i) where u, is the long wave velocity. Prior to wave 


breaking, bound long waves are forced by the groupiness of short waves in response to 
the changing short wave momentum. More sediment is put into suspension under the 
long wave trough stirred by the higher short waves than under the long wave crest where 
the short waves are lower. Since the velocity under the long wave trough is directed 
offshore, the resulting net sediment transport is offshore. Therefore, the short wave 
energy envelope and sediment transport are negatively correlated outside the surf zone 


(Figure 1). 





0 50 100 150 20 


Figure 1. | Depiction of wave groups normally incident on a beach with 
corresponding wave envelope (solid line), bound infragravity wave (dashed line) 
and free returning infragravity wave (dash-dotted line). Sediment is stirred under 

the short wave energy envelope. 


Inside the surf zone, the short waves break decoupling the forced long waves, 
which then become free waves. Inside the surf zone, the height of the short waves and 
momentum flux generally decrease due to wave breaking. The decreasing momentum 
flux is balanced by a hydrostatic pressure gradient of the mean water level resulting in set 
up of the mean water level inside the surf zone. The breaking wave height is controlled 
by the total water depth, which is slowly (compared to the time scale of the short waves) 
modulated by the infragravity waves. Higher short waves (increased stirring of 
sediments) are now associated with the crest of the long wave (onshore velocity), 
whereas now the sediment transport and infragravity waves are positively correlated 
moving sediment onshore inside the surf zone. 

How and to what degree infragravity waves affect bar generation is still in 
question. Thornton et al. (1996) and Gallagher et al. (1998) examined the beach at Duck, 
North Carolina for a variety of conditions and found that transport due to infragravity 
waves is generally subordinate in importance to asymmetry and the mean flow of 
longshore currents undertow. The infragravity sediment transport is not consistently in 
one direction, but alters as the cross-correlation changes throughout the surf zone. 
Ruessink (1998) and Roelvink and Stive (1989) all reason bar generation is the result of 
the varying degree of influence of multiple flow systems related to wave asymmetry, 


undertow, and infragravity waves. 
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Il. DESCRIPTION OF MODEL 


A. OVERVIEW 


Delft3D is a depth integrated one-, two-, or three-dimensional comprehensive 
modeling system. The model simulates flow, wave, sediments, and morphological 
progression. Initially Delft3D was developed for modeling 3-D flow in estuaries. This 
study will utilize a research version of the two-dimensional hydrodynamic and 
morphodynamic model. The research version of Delft3D has modified wave driver and 


sediment transport modules (Figures 2 and 3). 


Standard D3D Mode! Research Version D3D Mode! 











Figure 2. Delft3D Model Tree: The research model contains a restart file that takes 
data after 15 minutes and sends it to the initial flow to restart the computational 
loop. 
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Figure 3. Delft3D Standard and research model: The research model does all of the 
above in one module instead of 4 separate modules. All information from one 
process is available to all other processes after each time step. 


1. Wave Model 
Delft3D-WAVE is a model used to simulate the propagation and transformation 
of wave energy from given initial environmental conditions of waves and wind over 


arbitrary bottom depths. In the standard version, the HISWA model (Holthuijsen et al., 
8 


1989) is used to take into account wave generation due to wind, energy dissipation due to 
wave-bottom interactions, and nonlinear wave-wave interactions. In the research 
version, the waves are phase-averaged over high frequency swell, but the infragravity 
band waves are phase resolved. HISWA solves for cg and 8 using the initial conditions 
and bathymetry. The short wave energy, E,,, is solved through the energy flux balances 


given by: 


OE, e 0E,,c , cos ) 0E,,c, cos@ ) oa (1) 
ot ox oy 





where cg is the group velocity, 8 is the incidence angle with respect to the x-axis, x is the 
distance in the cross-shore, y is the distance in the alongshore, and D,, is the wave energy 
dissipation. The dissipation, or decrease in organized wave energy, is used as the input 
for the roller energy flux model described below. Dissipation is modeled as 


D,, =P,D, where P, is the probability that a wave is breaking and Dj is the dissipation 


rate in a breaking wave (Roelvink 1993). Both P, and D, vary on the time-scale of the 
wave groups. In the research model, the energy term of grouped short waves is used to 
calculate the forcing of the infragravity wave. The probability that a wave is breaking is 


a function of the energy of the waves, F, and total water depth, h, defined as: 





n/2 
E 
P(E,h )=1- — 2 
ee ~ fa - 


1 : : : are : 
where E,,, = g Pah >| p is the mass density of water, g is gravitational acceleration, and 


y is the wave breaking parameter, and n is a dissipation parameter corresponding to the 
randomness of the incident waves. The dissipation of a breaking wave is analogous to 
the dissipation of a bore (Battjes and Janssen 1978), D, = 20f,£ , which is a function of 
the energy of breaking waves, the frequency, /f,, at the peak of the wave spectrum and a 


calibration coefficient, a , set equal to one. The total dissipation is found by multiplying 


P, and D, ; 





D,=P,D,= |1-—exp {Fe 2of,,E (3) 


The breaking waves are assumed to form propagating bores, or rollers, inside the surf 
zone. Time is needed to convert organized kinetic energy and potential energy into 
small-scale, dissipative turbulent motion that results in a temporal lag between initiation 
of breaking and eventual dissipation of wave energy. The breaking wave energy is used 


as the input, or source, to the roller energy flux equation, 


OE, , O2E.C, cos(@ ) _ O2E cos(6) eee (4) 
ot ox oy 


where E’, is the kinetic energy of the roller, c, is the phase speed, and D, is the turbulent 


energy dissipation in the roller given by Nairn et al. (1990), 


_ 2gsin( B)E, 
ee a (5) 


where f is the angle of the roller face and is set equal to 0.1. Therefore, the only free 


parameter to be determined in the Delft3D-Wave driver is the wave breaking parameter 


(Cy: 


2. Flow Model 


Delft3D-FLOW calculates nonsteady flow forced by waves, tide, and wind. This 
module provides the hydrodynamic basis for morphological computations. The 
infragravity waves are solved in the time domain using the nonlinear shallow water 
equations to phase resolve bound and free, trapped (edge waves) and leaky infragravity 
waves. The flow module solves an unsteady shallow water equation for an 


incompressible fluid to obtain a two-dimensional (depth averaged) simulation. The 


10 


velocity field is composed of the Eularian velocity, uv” , and includes the short wave drift 
velocity (Stokes 1847), 


u=u" +u*cos@) (6) 





The Stokes drift, u* = “ , (Phillips, 1977) must be subtracted from the computed 
p 


velocities, u and v, to solve for transport velocities. M, includes the contribution of the 


roller in the mass flux. 


pss (E. + 2E, )cos@) (7) 


x 
Cc 


In this approach, the vertical momentum equation is reduced to the hydrostatic 
approximation, i.e. vertical accelerations are assumed to be small compared to the 
gravitational acceleration and thus neglected. The depth averaged continuity equation is 
given by: 

hu)  dalhv 


cee + —— = 


ot ox oy 





where u and v are the wave group velocities in the x and y directions. Also h=d+n, 
where / is the total water depth, 7 is the mean water level, and d is still water level. The 


momentum equations in the x and y- direction are: 








2 2; 
BI OO ON Fn OO Wh 0G ) 
Ot ox oy Ox ph ph ox dy 


F a 8? 
Bg OV, 5 OV 2. ONE, On Fy. 2520 WO 2s (10) 


a? a. Oy oe Ph OR. ee Oe 





where T, , are the x and y component of the bed shear stress, which varies at the time 


11 


scale of the wave groups, and v is the turbulent eddy viscosity. Fx and Fy are the x and y 


components of the wave- induced forces: 


aS, OS, ds, as, 
Fx= —+— Fy=|—"+—> (11) 
ox dy 








where S,, are the components of wave momentum, or radiation stresses. The terms of the 


momentum equations in the x- and y- direction (going from left to right) are local change 
in velocity, advection terms, pressure gradient, bottom stress, external forces (wind and 
waves), and turbulent mixing. 

Waves and currents are two important hydrodynamic factors that dominate the 
coastal zone. The model takes into account the generation of currents by waves, such as 
undertow and longshore currents. Using time dependent depth averaged shallow water 
equations allows shear instabilities of mean sheared currents to develop. 

The bed shear stress of currents is enhanced by waves in the model. Bed shear 
stress is described through a nonlinear wave-current interaction within the bottom 
boundary layer, where turbulent shear stress is proportional to the velocity squared, 
tou’ (Soulsby et al. 1993). Bed shear stresses for currents and waves alone are 


separately calculated and then combined. The current bed stress is found using the 


Manning formulation: Tt, = pC, lv ‘| where IU | is the magnitude of the depth averaged 


2; 
horizontal velocity, the drag coefficient C, = cold , and n is the Manning coefficient. The 


vn 
‘ ; 1 = - 
magnitude of the bed stress due to waves is computed by: T,, = 5 of,U,,, where U,, is 
the amplitude of the near-bottom wave orbital velocity and f |, is the wave friction factor. 


The Soulsby parameterization of the Fredsoe (1984) model combines wave-current flow 


and time-mean bed shear stress, t, = y(t, +t,,), where y =x{1+bx” (i-x)'| and 


Cc 


x =——*—~_ and J, p, and q are fitting coefficients. 
(.+t,,) 
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3. Sediment Transport Model 


The research model computes sediment transport on the same time scale as the 
flow. An advection/diffusion equation model is used to for sediment transport (Galapatti, 
1983). 

p) p) p) hC,, ~hC 


—hC+—hCu*® +—hCv*? = —t+ —_ (12) 
ot ox oy T 


where C is the sediment concentration, w, is sediment fall velocity, and 7, is the 


adaptation time for the diffusion of the sediment given by, 


T, = 0.05 (13) 
WwW 


The equilibrium sediment concentration is obtained by the Soulsby-van Rijn sediment 
transport formulation (Soulsby 1997). 


1 24 


2 \2 
Exe p(A,,.A,,) (a* 4+pF i + O18U ins —Uu (1 — 3.5m) (14) 
eq i e C cr 
d 


where A,,and A,, are the bed load coefficients which are a function of the sediment grain 
size, relative density of the sediment, and the local water depth (Soulsby 1997). 
Presently A,,and A, do not include additional stirring of sediment by infragravity 
motions (van Ryn, 1993). The infragravity motions are assumed to be implicit in the 


near-bed orbital motion. To include infragravity velocities in the sediment stirring 


A0 


te =) 
h 


C,,, 18 due only to the mean current in which z, is set to .006. mm is the local bed-slope. 


requires a recalibration of A,, and A. C, = where the drag coefficient, 


u,, 18 the critical threshold that the mean and orbital velocities must surpass to stir 


CF 


sediment. 7” and 7” are the mean (averaged over many wave groups) Eulerian 
13 


velocities that stir the sediment, and u,,,. is the combined wave breaking induced 


turbulence motion and near bed short wave velocity, 


Urns = | Use +0.5k, (15) 


The near bed, wave breaking induced turbulence, k,, is given by (Roelvink and Stive 


1989). 





k, = (16) 





rms,s 


where Hrs, hi is the short wave, root mean square wave height. 

The research model calculates asymmetry by applying stream function theory 
using the local wave height, depth, and period as input to the time- integrated contribution 
of the wave normlinearity. Due to the sediment transport associated with asymmetric 
short waves and the phase lags between sediment transport and wave group motions, 


transport has been averaged over the peak wave period, T.,, giving the wave asymmetry 


related parameterized velocity i , 


fo C., Ou; at 


— (17) 
Jo’ C., Oat 


where C/, is the equilibrium sediment concentration and uw, is the instantaneous near bed 
velocity (prime denotes intrawave variables). The wave asymmetry calibration factor, 


Q,, is obtained by optimizing computed and observed measurements. The wave 


w? 


asymmetry related velocity # is added to u”, the instantaneous Eulerian velocity. 


Bottom changes are obtained from continuity. 


14 


nal ee ee (18) 





where p, is the sediment density, 1, is the porosity, and T 


morp 


, 18 a morphological time 


step. 
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IV. DUCK94 EXPERIMENT 


The data used for this study were acquired during the Duck94 experiment from 
September | until October 31, 1994 at the U. S. Army Corps of Engineers Field Research 
Facility (FRF) on a barrier island in the Atlantic Ocean. Two bars having a well- 
developed inner bar during this time period characterize the site. Waves and currents 
were measured with a cross-shore array of 15 bottom mounted bidirectional 
electromagnetic current meters and pressure sensors extending from near the shoreline to 


4.5 meters depth (Figure 4). 


Sensor Positions 





“400 150 200 250 300 350 400 450 500 550 


Meters 


Figure 4. Sensor positions and bottom profile (relative to mean sea level). Location 
and elevation of numbered pressure and current sensors indicated by circles. 


Offshore directional wave spectra, giving wave height, direction, period, and water depth, 
were measured using a linear array of 9 bottom mounted pressure sensors in 8m depth 
about 900m offshore. Data were sampled at 2 Hz. Bathymetric surveys were conducted 
regularly throughout the experiment by measuring the elevation and position using an 
amphibious vehicle (CRAB). The area surveyed, referred to as the minigrid, was 550m 


alongshore by 420m offshore composed of 18 cross-shore profile lines that where 25 
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meters apart near the instruments and 50 meters apart elsewhere. Details of the model 


calibration are described in the model calibration section that follows. 
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V. MODEL CALIBRATION 


The accuracy of the modeling is determined by the model parameters chosen. 
The optimal breaking wave free parameter y is first determined for Deflt3D by 
examining a period of alongshore beach uniformity. A second time period of alongshore 
nor-uniformity is chosen to study the 2-D affects on hydrodynamics and by model inputs, 
such as wave angle and wind data. The roller dissipation and bottom friction for both 1- 
D and 2-D are examined and set in accordance with the values found by Reniers et al. (in 
preparation 2002). The focus then turns to beach morphology by comparing the model 


and observed morphology with different asymmetry coefficients, ©, . 


A. 1-D HYDRODYNAMIC CALIBRATION 

The wave part of the hydrodynamic model was calibrated with data chosen from 
September 16" through the 28'" (yeardays 259 — 271), when waves arrived at near 
normal incidence, wave refraction was small and did not effect dissipation. Directional 
wave data having three-hour averages were obtained from the FRF 8m linear array. Hig 
ranged between 0.2 and 1.8 m, peak periods ranged between 3.2 and 15.6 s, and the wave 
angle ranged between 50° from the North and —38° from the South with a mean wave 
angle, 8, slightly from the South at -3.9 °degrees (Figure 5). 

The only free parameter in the Delft3D-Wave driver is the wave breaking 
parameter y (eq.3), where nd=10, the manning coefficient n = 0.02, and o (eq.3) and 
B (eq.5) are set to 1 and .1 respectively. Adjusting 8 affects the dissipation of roller 
wave motion. By decreasing B , the model roller dissipation is decreased. The system for 
organized wave motion (wave height) is governed by y. If y is decreased, dissipation 


from wave breaking increases and H,ms approaches Himax. 
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Figure 5. September 16-28 (yeardays 259-271) Wave Conditions: 
Figure 6. (A) Significant wave height measured at 8m water depth. (B) Wave 
period at the peak of the spectrum. (C) Wave direction measured with waves 
propagating straight on shore having zero degrees, positive from the North and 
negative from the South. 


Optimal y values were determined by comparing measured and computed wave 


heights for 15 model runs varying the wave breaking parameter (y ) between .03 and .06. 


A linear regression was completed to minimize the error of the breaker parameter ? as a 


function of Hms using model skill as a criterion. The skill is a scale where one is the 


highest value and indicates there is no difference between the model calculated data and 


the measured data. 
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where subscripts m and c denote measured and calculated data, and i is the instrument 
number. 
The skill for y values from .3 to .6 were compared to H,ms, tide, and wave period. 


Y was found to increase as a function of Hms (Figure 6). The linear regression for the 


wave breaking parameter is y =.074H,+.3180. As Hms increased, the optimal y for 


the model also increased. There was a broad area of high skill were waves varied 
between .6 and 1.2 meters. But as the waves increased in height the range of high skill 
decreased (Figure 6). Overall the model is not sensitive to reasonable variation of the 


wave breaking parameter y . 
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Figure 7. Skill plot (Hrms vs. y ) 


1-D hydrodynamics results found model and observed Hrms agree with an 


average skill of .89 for all 15 sensors for a y value of .425 during a thirteen-day period 
(Figure 7). A constant value of .425 was chosen for y throughout because the value 
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maintained the highest skill over the broadest range of wave height. The H,,, determined 
from pressure sensors from near shore (p5) to off shore (p19) are shown in Figure 7. The 
model was less accurate during the largest waves when the significant wave height 


reached about 1.8m. As the skill plot for y shows, the larger wave heights have a higher 
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Figure 8. | Cross-shore model and observed Hrms where modeled H,ms (line) and 
observed H,ms (dotted) are compared at different pressure sensors in the cross- 
shore using Y = .425 and giving a mean skill of .89. 


B. 2-D HYDRODYNAMIC CALIBRATION 

The 2D flow field is calibrated using data from October 10‘ through the 21™ 
(yeardays 283 - 294) because this time period is strongly 2-D. The Hig ranged between 
0.4 and 2.1 m at the FRF linear array for three-hour averages. The mean period was 


approximated as the period of the first moment, T.,,,. T,,,,, was computed by integrating 


the spectral density, S, and frequency, f, over frequency and direction. 
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Mean periods ranged between 6 and 11.6 s with a mean of 8.5 s. The wave angle ranged 
from 21° from the North to —18.2° from the South having average mean wave angle, 8, 


slightly from the North at 1.2 °. An effective mean wave angle, 6 , was computed that 
gives the same radiation stress calculated using linear theory as that determined from the 


full directional spectrum E(f,8) (see Thornton and Guza, 1986). 


sin cos0- (21) 


© [> | 


C5... 7 
S,, = /JEC(f,0)—sin 8 cosOdf@ = 
. fe Cc r : 
Alongshore winds were acquired with an anemometer located 19.5 meters above 
mean sea level at the end of the FRF pier. The alongshore winds ranged from -5.16 to 
14.15 m/s with a mean velocity from the North at 4.8 m/s over the time period specified 
(Figure 8). 
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Figure 9. October 10-21 (yeardays 283-294) wave conditions. (A) Significant wave 
height measured at 8m water depth. (B) Wave period. (C) Wave direction 
measured with wave propagating straight on shore having zero degrees, positive 
from the North and Negative from the South. (D) Alongshore wind where 
positive is from the North and negative is from the South. 
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Daily minigrid bathymetry and monthly largegrid bathymetry were available and 
necessary to model the longshore current. The minigrid data from October 10, 11, 12, 14, 
20, and 21 were superposed onto the September large-grid data to better resolve 
bathymetry at instrument locations. An extended flow field was needed to prevent 
boundary generated disturbances from affecting flow in the area of interest. 

October 10-15 was a period of offshore bar motion during relatively high-energy 
conditions. The largest observed Hms profile changes during the Duck94 experiment 
occurred between October 10" and 20". Initially the beach was alongshore uniform as 
seen by the 10 October minigrid bathymetry. The 18 October survey clearly shows 
alongshore inhomogeneitics, which developed during the storm (Gallagher et al., 1998). 


Also a large rip channel developed in the minigrid during this time period (Figure 9). 
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Figure 10. Grid for 10 and 20 October. The minigrid area shows the large rip 
channel (y=1000m), nonuniformity of the alongshore bathymetry, and the pier 
area (y=500m). 


In the model, the largegrid was expanded to 2,500 meters alongshore maintaining 
auniform AY and to 930 meters cross-shore with a systematically increasing AX from 5 
meters starting at the shore to 30 meters offshore. The grid-size for the model was 
chosen considering accuracy and computational time. Originally the computational grid 
was with AX increasing offshore having 103 grid points and AY equal to 10 meters 
giving 179 grid points for a total of 18,437 computational grid points. By changing AX , 
while still increasing with distance offshore, to only 66 grid points and AY to 12.5 
meters with 147 grid points, the total grid points were decreased to 9,702 decreasing 
computational time from 53 hours to 24 hours of computer time for 24 hours of real time. 
Comparing the skill for wave height and currents during October 10", there was 


essentially no decrease in skill by changing the grid (Figure 10). 
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Figure 11. | Skill comparison using different grids: H,ms, u, and v model skill values 
for October 10th compared using a fine grid and a normal grid. 

Delft3D uses an implicit finite difference numerical scheme (Stelling, 1984) to 
solve the momentum, continuity, and transport equations on a staggered grid. The 
Alternating Direction Implicit (ADI) technique is used and each time step is split into two 
stages. This scheme implicitly solves the water levels and velocities in the x-direction in 
the first half-step and the ydirection terms in the second half-step. Since the model is 
working at a wave group structure scale, the model allows for larger spatial scales with 
similar accuracy as smaller space scales. 

The 2-D shallow water equations are time integrated using a Crank-Nicholson 
implicit finite difference approximation. The implicit numerical scheme makes the 
computations more stable, allowing for larger spatial steps while maintaining accuracy. 

When calibrating the 2-D hydrodynamic part of Delft3D, the wind contribution 
was examined. The model improved when wind is included as a boundary condition as 


an example. The October 10" model run was completed with wind and without wind 
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input for comparison. The mean wind on October 10" was 9.8 m/s consistently coming 
from the North. The model run with wind data shows an increase in mean skill for the 
entire day (Figure 11). The mean skill for longshore current with no wind was .60, while 


the mean skill for longshore current with the wind was 0.65 showing an improvement of 
8.3%. 
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Figure 12. Longshore current skill compared with and without meteorology data. 


Bottom roller dissipation and friction were examined and found consistent with 
Reniers et al. (in preparation 2002) using a roller dissipation coefficient, B = .1 and a 
Manning number of 0.02. If B is decreased, roller dissipation decreases and longshore 
current shifts toward the shore. The Manning number affects bed shear stress. If the 


Manning number is increased, friction increases, decreasing the model flow (Figure 12). 
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Figure 13. Longshore velocity and manning number changes. 


Significant improvement in the modeled longshore current occurred when 
boundary conditions were modified. Initially, the North and South boundaries of the grid 
were reflective (v=0). When the longshore current approached the boundary, it was 
deflected offshore, but instead of fading away as depth increased it recirculated in the 
computational domain. The circulation caused the alongshore current furthest from the 
shore to be in the opposite direction as the current near the beach, which was in the same 
direction as the wind and wave angle. Therefore the reflective boundary conditions were 
replaced with water level boundary conditions that allow the flow to pass. 

The overall longshore current velocity skill for the 10'" through the 21" of 
October for all sensors was -0.2 (Figure 13). The beginning of the period showed 
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Figure 14.__ Longshore modeled and observed velocities where modeled longshore 
current (line) and observed longshore current (dotted) are compared at different 
velocity sensors in the cross-shore giving a mean skill of -0.2 


reasonably high skill, falling off around October 14'" (yearday 287) and then increasing 
again to October 21" (yearday 294). The bathymetry during this period changed 
significantly as the storm moved in, going from alongshore uniform to non uniform with 
a large rip channel. The model performed better before the storm. Also due to the storm, 
the only bathymetry available was October 10, 11, 12, 14, 20, and 21. The model 
performed better when daily bathymetry was available (Figure 14). The longshore model 
velocity skill was affected more by the lack of bathymetry than the cross-shore velocity. 
The cross-shore skill improved when the rip current dominated the flow (yearday 288) 


due to the sensitivity of the skill factor when measurements are close to zero. 
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Figure 15. | Cross-shore and longshore skill for 10-21 October with negative skill 
values are not shown. 


Examining 19-20 October (yeardays 292-293) for sensors in the trough and on the 
bar (v5, v22, and v13) the modeled longshore velocities are greater than the near- zero 
observed longshore velocities. The wave angle during this time was near zero. When the 
bathymetry was alongshore uniform, small wave angle (slightly varying from shore 
normal) did not affect alongshore and cross-shore currents. However, small wave angles 
resulted in strong variability in alongshore currents during times of complex bathymetry, 


such as after the storm. 


To show the sensitivity of the model, the 20" of October was examined because 
acceptable 2-D bathymetry information was available, the wave angle was shore normal, 
and the alongshore wind was minimal. Adjusting the wave angle by adding 2.5°, which 
is about the measurement error, the model was compared with a normal wave angle to 
determine how the computed results varied. As expected, H;ms did not change 


appreciably. On the other hand, the cross-shore current skill increased while the 
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longshore current skill decreased. The significant skill changes in opposite direction with 
such a small change in wave angle shows that currents and circulations in a rip channel 


are sensitive to small changes in wave conditions (Figure 15). 
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Figure 16. Skill comparison with 2.5° increased wave angle. 


Looking closely at the end of the period, the skill improved from the 20" to the 
21“ of October (Figure 13), but at the same time, the bathymetry and alongshore wind 
were almost identical. The question still remains why there are these differences between 


the two days. 


Another consideration for skill variations at this test location is the near proximity 
of the FRF 800m pier. The pier scour hole and pilings can influence the waves and 
currents and resulting morphology within the minigrid. The model does not consider the 


turbulent eddies from the pier or their effects on circulation. 
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VI. 2-DMORPHOLOGY ASSESSMENT 


The changes in morphology were examined for October 12'" through the 16" by 
comparing with the minigrid bathymetry. Asymmetry as determined by the ,, value, is 
important in determining net sediment transport. O,, was examined by varying its values 
between .1 and .4. If @, is decreased, onshore sediment transport is decreased, moving 


the bar offshore as sand erodes. The asymmetry is also a function of the wave period, 
where as the wave period increases, the asymmetry increases. The wave period increased 


throughout the period considered. 


Examining a cross-shore profile transect at 180m North of the velocity sensors, 


increasing O, significantly moved sediment onshore over the five day period (Figure 16). 


Cross shore transect y=1100m after 120 hours 
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Figure 17. Cross-shore sediment transport while increasing a, at y = 1100m. 
However, in the rip channel where the velocity sensors are located, increasing Q,, 


actually increased erosion and the sandbar moved further offshore and widened (Figure 


17). The model deposited sand further offshore than in actuality. As asymmetry was 
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increased, more sand moved onshore, increasing the set-up and pressure gradients on 
either side of the channel. This increased forcing, increased the magnitude of the rip 
current. The model does not include wave-current interaction throughout the flow field, 
except for the bottom shear stress, so as the rip increased in strength there was no 
opposing force to slow it down. In reality incoming waves opposing the outgoing rip 
current cause the sand to be deposited closer to the shore than what the model predicts. 


The rip current channel was widened and deepened as ©, increased. 
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Figure 18. Cross-shore sediment transport while increasing o,, at y = 930m. 
The modeled bathymetry for all @, values eroded more alongshore near the 


beach than the observed bathymetry. This indicates that the alongshore current was 


stronger right near the shore than observed. By comparing the modeled and measured 
minigrids after five days, it appears @,,=.25 most closely simulates the 2-D changes in 


bathymetry (Figure 18). 


34 


Initial Minigrid Final Minigria 
1100 5 





1050 
1000 
950 
900 
850 


800 








750 1 oe, 
400 200 300 400 400 200 300 400 


1100 
1050 
1000 
950 
900 
850 


800 





750 Sn S 
100 200 300 400 


1100 
1050 
1000 
950 
900 
850 
800 











XN 750 % 
100 200 300 400 100 200 300 400 





Figure 19. | 2-d bathymetry showing rip channel. 
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The effects of the scour hole created by the FRF pier was examined next. A 
smoothed bathymetry was generated removing the pier scour hole to examine the 
circulation and morphology. It appears that the model leaves a large deposit of sediment 
at the Southern boundary of the minigrid when the pier is included. Also the rip channel 
has an orientation from the Southwest to the Northeast with the pier scour hole included 
in the model. Removing the effects of the pier, the orientation of the channel changes to 
the Northwest to the Southeast with a larger sediment deposit on the Northern boundary 
of the minigrid. From this it is clear the pier plays a significant role in the sediment 


dynamics in the area. 
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VU. CONCLUSIONS 


Nearshore morphology and currents were measured nearly continuously for two 
months near Duck, North Carolina where September 16-28 and October 10-21 were two 
time periods identified to model. Wave height time series outside the surf zone show 
considerable variability in wave height associated with wave generation by local winds 
and distant storms. Strong winds from Nor’esters build waves increasing the radiation 
stresses and turbulent stirring of sediment. Due to near shore circulations, a rip channel 
developed and deepened due to the increased offshore forcing. Nearshore modeled H-ms, 
currents, and morphology were compared to data by the computed Delft3D 2-DH model 
utilizing shallow water equations to phase resolve the mean infragravity motions in 


combination with an advection diffusion equation for the sediment transport. 


Energy dissipation due to wave breaking is extremely complex. The wave 
breaking parameter, y , utilized by Roelvink (1993), is critical in the decay of wave 
groups. Gamma is found to be a function of wave height having a concentration of high 
skill around gamma=.425 for a broad range of wave heights. Interestingly, Roelvink 
(1993) found a higher optimal value of gamma using lab data of .57. From the 
verification presented (Figure 6), we may draw the conclusion that the model is relatively 


insensitive to varying gamma values. Concluding, y does not need to be adjusted to fit 


field data. 


The current bed shear stress is related to the horizontal velocity scaled by the drag 
coefficient. The drag coefficient is proportional to the Manning number squared. The 
tests indicate the model predictions vary with variations of the Manning number but are 
not overly sensitive. The model was sensitive to small changes in wave angle when 
waves were near shore normal and the bathymetry was complex. Significant changes in 
skill took place with small adjustments to the wave angle. As it is difficult to accurately 


measure wave angle this is an important circumstance in the model assessment. 


The asymmetry coefficient primarily responsible for moving sediment onshore 


due to cross-shore wave velocity asymmetry strongly affects model sediment transport. 


ea 


Qualitatively there is reasonably good agreement between measured and predicted 
sediment distribution over the 2D minigrid area whena®,, = .25 or @ = .3. The rip 
channel was modeled in roughly the same location and with a clearly developed channel. 
Further study is required to calibrate the model. 

Man made objects affect the wave, flow, and sediment dynamics. The model 
does not resolve sub grid objects like the pier pilings. Man made objects can produce 
turbulent eddies which are on a time and special scale smaller than the model can resolve. 
A possibility would be to explore is the idea of parameterizing the bottom stress of the 


model to take into account the turbulence by a man made object. 
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